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Abstract 

A  full-spectral  third-generation  ocean  wind-wave  model  (Wavewatch-III)  implemented  in  the  South  China  Sea  is  used  to 
investigate  the  effects  of  the  wave  boundary  layer  on  the  drag  coefficient  and  the  sea-to-air  transfer  velocity  of  dimethylsulfide 
(DMS)  during  passage  of  Typhoon  Wukong  (September  5-11,  2000)  with  a  maximum  sustained  wind  speed  of  38  m  s-1.  The 
model  is  driven  by  the  reanalyzed  surface  winds  (1°  x  1°,  four  times  daily)  from  the  National  Centers  for  Environmental  Prediction. 
It  is  found  that  the  wave  boundary  layer  evidently  enhances  (16.5%)  the  drag  coefficient  (in  turn  increases  the  momentum  flux 
across  the  air-sea  interface),  and  reduces  (13.1%)  the  sea-to-air  DMS  transfer  velocity  (in  turn  decreases  the  sea-to-air  DMS  flux). 
This  indicates  the  possibility  of  important  roles  of  wave  boundary  layer  in  atmospheric  DMS  contents  and  global  climate  system. 
©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  dominant  natural  source  of  sulfur  to  the  atmo¬ 
sphere  is  the  oceanic  dimethylsulfide  (DMS)  (Bates 
et  al.,  1992;  Gondwe  et  al.,  2003).  In  the  atmosphere, 
DMS  is  oxidized  to  sulfuric  and  methanesulfonic  acids 
which  condense  to  form  new  aerosol  particles  and/or  to 
add  mass  to  existing  particles.  These  particles  can  affect 
the  Earth’s  radiation  budget  by  scattering  solar  radiation 
back  to  space  and  altering  the  properties  and  lifetimes  of 
clouds.  In  estimation  of  the  Earth’s  radiation  budget,  the 
sea-to-air  DMS  flux  must  be  included  in  the  chemical 
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transport  and  climate  models  in  order  to  accurately 
calculate  aerosol  radiative  forcing.  The  sea-to-air  DMS 
flux  depends  on  airside  and  waterside  transfer  velocities. 

A  moving  tropical  cyclone  is  an  intense  source  of 
surface  wind  stress  that  produces  strong  ocean  waves 
which  changes  the  surface  momentum  flux  and  sea-to- 
air  DMS  transfer  velocity.  Up  until  now,  effect  of  wave 
boundary  layer  (WBL)  on  the  surface  momentum  flux 
through  drag  coefficient  CD  and  roughness  length  z0  has 
been  investigated  over  small  areas  of  the  sea  or  in  wave 
tanks  (Hwang,  2005),  but  not  over  regional  seas  with 
high  waves  forced  by  strong  winds  such  as  by  tropical 
cyclones.  However,  effect  of  WBL  on  the  sea-to-air 
DMS  transfer  velocity  has  not  yet  been  investigated.  In 
this  paper,  we  study  the  effects  of  WBL  on  (CD,  z0)  and 
sea-to-air  DMS  transfer  velocity  in  the  South  China 
Sea  (SCS)  during  the  passage  of  Typhoon  Wukong 
(September  4-10,  2000). 
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A  fully  spectral  third-generation  ocean  wind-wave 
model,  Wavewatch-III  (henceforth  denoted  as 
WWATCH),  is  used  to  simulate  the  WBL  parameters. 
The  effects  of  WBL  can  be  easily  identified  by  the 
differences  in  (CD,  z0)  and  DMS  transfer  velocity  between 
with  and  without  WBL.  The  outline  of  this  paper  is  as 
follows:  A  description  of  SCS,  Typhoon  Wukong  2000, 
and  WWATCH  model  is  given  in  Sections  2,  3,  and  4. 
Drag  coefficient  and  roughness  length  are  given  in  Section 
5.  The  effects  of  WBL  on  the  drag  coefficient,  roughness 
length  and  DMS  transfer  velocity  are  presented  in  Sections 
6  and  7.  In  Section  8  the  conclusions  are  presented. 

2.  South  China  Sea 

The  SCS  is  a  semi-enclosed  tropical  sea  located 
between  the  Asian  land  mass  to  the  north  and  west,  the 
Philippine  Islands  to  the  east,  Borneo  to  the  southeast, 
and  Indonesia  to  the  south,  a  total  area  of  3. 5 TO6  km2. 
Its  southern  border  is  3°  S  between  South  Sumatra  and 
Kalimantan  (Karimata  Straits),  and  its  northern  border  is 
the  Strait  of  Taiwan  from  the  northern  tip  of  Taiwan 
to  the  Southeast  coast  of  China.  All  of  the  straits  are 
shallow  except  Luzon  Strait  whose  maximum  depth  is 
1800  m.  The  elliptical  shaped  central  deep  basin  is 
1900  km  along  its  major  axis  (northeast-southwest)  and 
approximately  1100  km  along  its  minor  axis,  and  ex¬ 
tends  to  over  4000  m  deep. 

The  SCS  can  be  viewed  as  a  distinctive  ecosystem 
due  to  its  boundaries  of  archipelagoes  and  peninsulas, 
dotted  by  small  islands  and  coral  reefs;  the  striking 


variation  in  its  sea  floor  characteristics,  averaging  1 00  m 
deep  on  the  continental  Sunda  shelf  and  over  5000  m  in 
the  Philippine  basin;  and  its  unusual  monsoon  weather 
patterns  of  reversing  summer  and  winter  rains  and 
winds.  The  northeast  monsoon  between  December  and 
February  and  the  southwest  monsoon  between  June  and 
August  change  the  surface  water  circulation  pattern. 
Geology  and  climate  combine  to  produce  a  remarkable 
amount  of  biological  diversity  and  immense  genetic 
resources  in  the  SCS.  Extensive  coral  reefs  support 
several  thousand  different  species  of  organisms  and  play 
an  important  part  in  buffering  wave  impact  on  beaches, 
thus  reducing  erosion.  The  highly  variable  wind  systems 
and  complicated  topography  make  SCS  a  perfect  loca¬ 
tion  for  the  investigation. 

3.  Typhoon  Wukong  2000 

Typhoon  Wukong  2000  formed  as  a  tropical  depres¬ 
sion  west  of  Luzon  Island  at  06UTC  on  4  September.  It 
remained  almost  stationary  until  12UTC  on  5  September 
and  then  made  an  anticlockwise  turn  to  the  west  (Fig.  1). 
During  the  turn  the  depression  was  developed  into  a 
tropical  storm  at  00UTC  on  6  September.  It  further 
intensified  to  attain  modest  typhoon  intensity  at  1 8UTC 
on  7  September  and  reached  its  peak  with  estimated 
maximum  sustained  winds  of  38  m  s-1  at  06UTC  on 
8  September  (Fig.  2).  Weakening  gradually,  Wukong 
skirted  around  the  southern  coasts  of  Hainan  Island  on  9 
September  and  made  landfall  on  the  northern  part  of 
Vietnam  at  around  04UTC  on  10  September.  After 
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Fig.  1.  Track  of  Typhoon  Wukong  2000  during  September  4-10,  2000. 
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Fig.  2.  GMS-5  visible  image  of  Tyhoon  Wukung  at  003 1Z  on  8  September  2000  when  the  cyclone  was  located  about  390  km  east  of  Hainan  Island 
and  tracking  westward.  Evident  outflow  is  seen  in  all  quadrants  as  well  as  a  clear  but  irregular  eye. 


the  landfall  it  was  weakened  into  a  tropical  depression  in 
the  northeastern  part  of  Thailand  at  12UTC  on  10 
September. 

4.  W WATCH  model 

4.1.  Description 

WWATCH  was  developed  for  regional  sea  wave 
prediction  at  the  Ocean  Modeling  Branch  of  the  Envi¬ 
ronmental  Modeling  Center  of  the  National  Centers  for 
Environmental  Prediction  (NCEP).  It  was  built  on  the 
base  of  Wavewatch-I  and  Wavewatch-II  as  developed 
at  the  Delft  University  of  Technology  and  NASA  God¬ 
dard  Space  Flight  Center,  respectively  (Tolman,  1999). 
WWATCH  implemented  in  the  SCS  has  been  evaluated 
using  the  significant  wave  height  (SWH)  data  collected 
by  satellites  with  altimetry,  such  as  TOPEX/POSEIDON 
(T/P)  (Chu  et  al.,  2004).  The  SWHs  from  WWATCH  are 
compared  to  that  from  T/P  at  the  satellite  crossover 
points  in  SCS.  The  model  errors  of  SWH  have  Gaussian- 
type  distribution  with  small  mean  value  of  0.02  m 
(almost  no  bias).  The  model  errors  are  comparable  to  the 
T/P  altimeter  accuracy  (0.5  m)  in  the  central  SCS  and 
smaller  than  the  T/P  altimeter  accuracy  in  the  northern 
and  southern  SCS,  which  indicates  the  capability  of 
WWATCH  for  the  SCS  wave  simulation. 


Let  (x,  y)  be  the  Cartesian  coordinates  with  (i,  j)  the 
unit  vectors.  The  wave  spectrum  F  is  generally  a  func¬ 
tion  of  space  (x=  xi+yj),  time  (f),  and  phase  parameters 
(cr,  k,  9,  co), 

F  =  F(k,9,<r,co;x,t),  (1) 

where  k  is  the  wave  number  (unit:  m  '  );  0  is  the  direction 
(unit:  degree);  cr  is  the  intrinsic  frequency  (unit:  s_  *);  and 
co  is  the  absolute  frequency  (unit:  s_l).  The  individual 
spectral  components  are  assumed  to  satisfy  the  linear 
wave  theory  (locally)  and  to  follow  the  dispersion  relation, 

o2  =  gk  tank  kd ,  co  =  cr  +  kU ,  (2) 

where  d  is  the  mean  water  depth;  U  is  the  (depth-  and 
time-averaged)  current  velocity.  When  U  is  given,  only 
two-phase  parameters  among  (cr,  k ,  9 ,  co)  are  independent. 
Many  wave  models  use  the  frequency-direction  (cr,  9)  as 
the  independent  phase  variables.  However,  WWATCH 
uses  the  wavenumber-direction  (k,  9)  as  the  independent 
phase  variables. 

Without  currents,  the  energy  of  a  wave  package  is 
conserved.  With  currents  the  energy  of  a  spectral  com¬ 
ponent  is  no  longer  conserved  (Longuet-Higgins  and 
Stewart,  1961),  but  the  wave  action  spectrum, 

N (k,  9;  x,  t)  =F(k ,  9;  x,  t)  / cr, 
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is  conserved  (Whitham,  1965;  Bretherthon  and  Garrett, 
1968).  In  WWATCH,  the  basic  equation  is  for  the  wave 
action  spectrum. 

4.2.  Model  setting 

WWATCH  has  two  types  (mandatory  and  optional) 
of  model  switches  for  users  to  choose.  Table  1  lists  the 
model  setting  and  optional  switches  for  this  study.  For 
example,  spatial  and  spectral  grids  are  user-defined;  the 
quickest  propagation  scheme  is  selected  with  the  dis¬ 
persion  correction  (Booij  and  Holthuijsen,  1987);  non¬ 
linear  interactions  are  included;  and  the  source  term 
parameterization  includes  wind  input,  nonlinear  wave- 
wave  interaction,  dissipation,  and  wave-bottom  interac¬ 
tion  (Tolman  and  Chalikov,  1996). 

4.3.  Discretization 

The  model  is  implemented  for  SCS  (0°  to  25°  N, 
105°- 122°  E)  with  realistic  bathymetry  data  from  the 
Naval  Oceanographic  Office  (i.e.,  DBDB5)  and  a  regular 
longitude-latitude  grid  (  Ax  =  Ay =25  km).  The  wave 


Table  1 

Model  setting  for  this  study 


Switch 

parameters 

Characteristics 

DUM 

Dummy  to  be  used  if  WWATCH  is  to  be  installed  on 
previously  untried  hardware 

LRB8 

8  byte  words 

SHRD 

Shared  memory  model,  no  message  passing 

SEED 

Seeding  of  high-frequency  energy 

GRD1 

Settings  directly  hardwired  to  user-defined  spatial 
grids  (spherical  coordinate  with  0.25°  grids) 

SP1 

User-defined  spectral  grids. 

PR2 

Ultimate  quickest  propagation  scheme  with  Booij  and 
Holthuijsen  (1987)  dispersion  correction 

ST2 

Tolman  and  Chalikov  (1996)  source  term  package 

STAB2 

Enable  stability  correction  for  Tolman  and  Chalikov 
(1996)  source  term  package 

NL1 

Nonlinear  interaction  (DIA) 

BT1 

JONS  WAP  bottom  friction  formulation 

WIND2 

Approximately  quadratic  interpolation 

CUR2 

Approximately  quadratic  interpolation 

ol 

Output  of  boundary  points  in  grid  preprocessor 

o2 

Output  of  the  grid  point  status  map  in  grid  preprocessor 

o2a 

Generation  of  land-sea  mask  file  mask.ww3  in  grid 
preprocessor 

o3 

Additional  output  in  loop  over  fields  in  field  preprocessor 

o4 

Print  plot  of  normalized  1  -D  energy  spectrum  in  initial 
conditions  program 

o5 

2-D  energy  spectrum 

06 

Spatial  distribution  of  wave  heights  (not  adapted  for 
distributed  memory) 

o7 

Echo  input  data  for  homogeneous  fields  in  generic  shell 

direction  ( 9 )  is  also  regularly  gridded  with  A0  =  15°.  The 
wavenumber  grid  spacing  is  determined  by  the  frequen¬ 
cy  intervals  (total  25) 

om+x=Xaom  m  =  0,1,. ..,24,  (3) 

with 

XG  =  1.1,  (j0  =  0.0418.  (4) 

Two  time  steps  (300  s  and  100  s)  are  used  in 
WWATCH  to  reach  computational  efficiency:  (a)  global 
time  step  (300  s)  for  the  propagation  of  the  entire 
solution,  (b)  spatial  time  step  (300  s)  for  the  spatial 
propagation,  (c)  spectral  time  step  (300  s)  for  intra- 
spectral  propagation,  and  (d)  source  time  step  (100  s)  for 
the  source  term  integration. 

4.4.  Model  integration 

WWATCH  is  integrated  from  the  JONSWAP-1973 
wave  spectra  (Hasselmann  et  al.,  1980)  with  four-time 
daily  reanalyzed  NCEP  winds  (1°  x  1°)  from  January  1, 
2000  to  31  December  2000.  The  output  of  WWATCH 
consists  of  the  traditional  frequency-direction  spectrum 
F(cr,0,  x,  t),  which  is  calculated  from  wavenumber- 
direction  spectrum  F(k ,  6 ,  x,  t )  using  the  Jacobean 
transformation.  The  peak  frequency  (&>p)  and  in  turn  the 
peak  phase  speed  (cp)  can  be  determined  for  (x,  t)  from 
the  spectrum  F(cr,0,  x,  t). 

4.5.  Model  verification 

SWH  is  chosen  for  model  verification  (Chu  et  al., 
2004).  The  model-data  comparison  is  conducted  at  all 
20  crossover  points  in  SCS.  Each  one  contains 
approximately  72  pairs  of  modeled  ( Hm )  and  observed 
( H0 )  SWH  data  in  2000.  The  total  number  of  pairs  is 
1330.  The  model  error  is  represented  by  the  difference 
between  the  modeled  and  observed  SWH  at  the 
crossover  points,  A H=Hm~H0. 

Over  the  whole  SCS,  WWATCH  has  very  low  bias 
(-0.01  to  0.04  m)  in  predicting  SWH  with  a  maximum 
(positive  bias)  value  of  0.04  m  in  March  and  a  minimum 
(negative  bias)  value  of  -0.01  m  in  April.  The  root- 
mean-square  error  has  a  minimum  value  of  0.39  m  in 
March  and  a  maximum  value  of  0.48  m  in  December. 
The  model  errors  are  comparable  to  the  T/P  altimeter 
accuracy  (0.5  m)  in  the  central  SCS  and  smaller  than  the 
T/P  altimeter  accuracy  in  the  northern  and  southern 
SCS,  which  indicates  the  capability  of  WWATCH  for 
SCS  wave  simulation.  Interested  readers  are  referred  to 
Chu  et  al.  (2004). 


126 


PC.  Chu,  K.-F.  Cheng  /  Journal  of  Marine  Systems  66  (2007)  122-129 


5.  Drag  coefficient  and  roughness  length 

The  mean  wind  profile  is  represented  by 

«(z)~£ln(T),  (5) 

K  \Z0J 

where  k  =  0.4  is  the  von  Karman  constant  and  u*  is  the 
air  friction  velocity,  which  is  defined  by 

CD  =  ul/u (6) 

where  CD  is  the  drag  coefficient  of  the  ocean  surface.  CD 
was  originally  considered  as  a  constant  (e.g.,  Kraus,  1 972). 
However,  as  the  quantity  and  quality  of  measurements 
have  improved,  it  became  evident  that  CD  tends  to  in¬ 
crease  with  increasing  wind  speed  (e.g.,  Garratt,  1977). 

Using  Eqs.  (5)  and  (6),  the  roughness  length  z0 
becomes 


z o  =  zr  exp 


—)■ 

VCdJ 


(?) 


the  atmospheric  boundary  layer  above  the  sea,  whose 
structure  is  influenced  directly  by  surface  waves.  Within 
the  WBL,  some  portion  of  momentum  transfer  results 
from  wave-produced  fluctuations  of  pressure,  velocity, 
and  stresses.  The  drag  coefficient  CD  at  the  reference 
height  zr  is  calculated  by  (Charlikov,  1995) 

CD  =  K1[R-lnCD}2,  (10) 

/?  =  In  ( — — — ^  ,  (11) 

\yyj^FrJ 


where  y  =  2  is  a  constant;  jup  is  the  conventional  non- 
dimensional  energy  level  at  peak  frequency  (&>p). 
Estimation  of  the  drag  coefficient  thus  requires  an 
estimate  of  the  energy  level  jup  at  peak  frequency,  which 
is  parameterized  by  (Janssen,  1989), 

Hp  =  0.57  (~)  ,  (12) 


The  roughness  length  is  also  expressed  in  the 
dimensionless  form 

z0*=z0g/ul,  (8) 

where  g  is  the  gravitational  acceleration,  and  the  dimen¬ 
sionless  roughness  z0 *,  is  generally  referred  to  as  the 
Chamock  parameter.  Without  WBL,  the  Chamock  pa¬ 
rameter  is  considered  as  a  constant  (e.g.,  Chamock,  1 955), 

z0*  =  0.0144.  (9) 

For  the  given  wind  speed  at  the  reference  height 
zr  (zr=  1 0  m  in  this  study),  the  drag  coefficient  CD,  rough¬ 
ness  length  z0,  and  the  friction  velocity  u *  can  be  obtained 
through  the  iteration  using  Eqs.  (6)-(9).  The  iteration 
stops  when  the  relative  change  of  the  friction  velocity  is 
smaller  than  a  prescribed  criterion  (10-3).  Such  iterations 
are  performed  during  the  model  initialization,  but  are  not 
necessary  during  the  actual  model  mn  as  w*  changes 
slowly  (Tolman  and  Chalikov,  1996;  Tolman,  1999).  The 
effect  of  the  atmospheric  instability  on  the  friction 
velocity  is  parameterized  using  an  effective  wind  speed 
we(Tolman  and  Booij,  1998),  which  depends  on  the 
surface  air  and  sea  temperature  difference. 

6.  Effect  of  WBL  on  drag  coefficient  and 
roughness  length 

Field  data  indicate  convincingly  that  CD  is  sea  state 
dependent  (Hwang,  2005).  The  WBL  is  the  lower  part  of 


Obviously,  jup  also  depends  on  the  drag  coefficient 
CD.  Thus,  the  iteration  should  be  conducted  to  get  (CD, 
z0,  uf)  using  Eqs.  (1 0)— (12).  During  the  iteration,  the 
peak  phase  speed  cp  in  Eq.  (12)  is  obtained  from  the 
WWATCH  simulated  wave  spectmm  (F). 

Fig.  3  shows  the  scatter  diagram  of  the  Chamock 
parameter  z0*  versus  the  nondimensional  peek  angular 
frequency  cop*  (=copu*/g)  with  WBL  during  the  passage 
of  Typhoon  Wukong  (4-10  September  2000).  The 
points  are  deviated  a  lot  from  the  horizontal  line, 


Fig.  3.  Comparison  of  the  Chamock  parameter  z0*  without  WBL  (solid 
line)  and  with  WBL  (dots)  for  the  northern  SCS  (15°-25°  N)  during 
5-10  September,  2000.  Here,  cOp^cOpU^/g  with  cop  the  angular  peek 
frequency. 
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zo*  =  0.0144  (without  WBL),  which  indicates  the  large 
effect  of  WBL  on  the  roughness  length.  The  horizontally 
averaged  drag  coefficient  (Fig.  4)  over  the  northern  SCS 
(15°-25°  N)  is  evidently  larger  with  WBL  (C’))v),  solid 
curve)  than  without  WBL  (Cd\  dashed  curve). 

The  relative  increase  of  the  drag  coefficient  due  to 
WBL  is  calculated  by 


6  = 


EE 

i  j 


E  E  do\u) 


(13) 


where  the  summation  is  over  the  total  number  of  grid 
points  in  the  northern  SCS  (15°-25°  N).  The  parameter 
8  generally  increases  with  the  intensity  of  Typhoon 
Wukong  (Fig.  5).  As  the  low  pressure  of  the  typhoon 
reaches  its  minimum  value  of  955  hPa  on  September  8, 
8  has  its  maximum  value  of  16.5%.  A  recent  study 
shows  reduced  drag  coefficient  for  high  wind  speeds 
(>40  m  s-1)  in  tropical  cyclones  (Powell  et  al.,  2003). 
However,  Wukong  2000  is  a  modest  typhoon  with 
maximum  sustained  winds  of  38  m  s_1.  The  simulated 
enhancing  WBL  effect  (on  drag  coefficient)  with  wind 
speed  during  the  passage  of  Wukong  2000  does  not 
conflict  with  Powell  et  al.’s  (2003)  results. 


7.  Effect  of  WBL  on  DMS  transfer  velocity 

The  sea-to-air  DMS  flux  exhibits  more  complexity 
than  the  bulk  aerodynamic  relations  described  in  the 
previous  sections.  Turbulent  transfer  in  the  atmosphere 
above  the  surface  plays  a  key  role,  but  chemical  proper¬ 
ties  of  the  surface  and/or  reactions  with  sea  spray  are 
also  important.  For  chemical  compounds  which  react 


Fig.  5.  Temporally  varying  relative  increase  of  CD  between  with  and 
without  WBL  (solid  curve)  for  the  northern  SCS  (15°-25°  N)  and 
typhoon  intensity  (dashed  curve)  represented  by  minimum  pressure. 
Note  that  the  increase  can  reach  16.5%  on  8  September  2000  when  the 
typhoon  Wukong  is  strongest. 


with  sea  spray  aerosols  (e.g.,  ammonia  and  nitric  acid), 
the  constant  flux  layer  assumption  may  not  be  valid;  and 
additional  equations  must  be  considered.  For  chemical 
compounds  with  reaction  time  scales  longer  than  the 
surface  layer  turbulence  time  scales,  the  constant  flux 
layer  assumption  may  be  invoked.  This  applies  to  gases 
such  as  DMS.  It  is  noted  that  the  photo  reactions  (same 
time  scale  as  the  turbulence)  are  not  considered  here. 
The  sea-to-air  DMS  flux  ( H)  is  characterized  by  transfer 
velocities  k a  and  kw  (McGillis  et  al.,  2000), 


H  = 


1  H-  otk-w  j  k% 


(Cw  oi(2a), 


(14) 


where  Ca  and  Cw  are  the  bulk  concentrations  in  the  air  and 
water;  a  is  the  Ostwald  solubility  coefficient  for  DMS, 


a  =  exp 


'3525 

~Y~ 


-9.464 


(15) 


Fig.  4.  Comparison  of  temporally  varying  horizontal  averaged  CD 
without  WBL  (dotted  curve)  and  with  WBL  (solid  curve)  for  the 
northern  SCS  (15°-25°  N)  during  5-10  September,  2000. 


where  Tis  temperature  (in  K).  Waterside  transfer  velocity 
( kw )  is  parameterized  by  (Jahne  et  al.,  1987;  McGillis 
et  al.,  2000) 

kw  =  \J  FfPw  /T1  Scnu* ,  (16) 

where  (pa=  1.29  kg  m-3,  pw=  1025  kg  m-3)  are  the 
characteristic  air  and  water  densities;  P  is  the  dimension¬ 
less  Stanton  number  showing  the  resistance  to  momentum 
transfer  across  the  water-side  viscous  boundary  layer;  and 

Sc  =  v/D ,  (17) 

is  the  Schmidt  number.  Here,  v  is  the  kinematic  viscosity 
for  momentum,  and  D  is  the  diffusion  coefficient  for  gas 
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transfer.  The  diffusion  coefficient  of  DMS  is  calculated  by 
(Saltzman  et  al.,  1993) 


D  =  1.1  x  10  2  exp 


1896 


(unit :  cm2s  1). 


(18) 


The  Schmidt  number  for  DMS  at  300  K  is  taken  as 
(Saltzman  et  al.,  1993), 

Sc  =  120.  (19) 

The  kinetic  viscosity  of  DMS  can  be  calculated  from 
D  and  Sc  using  Eqs.  (17)— (19). 

Airside  transfer  velocity  (ka)  is  calculated  by 

K  =  kn2o(M /My[2o)  ^2,  (20) 

where  M(=  62. 129)  and MHi 0  (=18.015)  are  the  molec¬ 
ular  weights  for  DMS  and  H20;  ku2 o  is  the  water  vapor 
transfer  velocity  which  is  calculated  by 

^h2o  =  659ur.  (21) 

In  fully  rough  regime  of  the  sea  state  such  as  the  case 
discussed  in  this  paper  (wave  generation  during  the 
typhoon  passage),  the  exponent  in  Eq.  (16)  is  given  by 
/2  =  0.58;  and  the  Stanton  number  is  represented  by 
(Jahne  et  al.,  1987) 

/?  =  0.55tfey4,  (22) 

where 


Fig.  7.  Temporally  varying  relative  decrease  of  waterside  transfer 
velocity  kw  between  with  and  without  WBL  (solid  curve)  for  the 
northern  SCS  (15°-25°  N)  and  typhoon  intensity  (dashed  curve) 
represented  by  minimum  pressure.  Note  that  the  relative  decrease  of 
waterside  transfer  velocity  follows  the  intensity  of  Typhoon  Wukong. 
It  reaches  - 1 3 . 1  %  on  8  September  2000  (strongest  strength  of  typhoon 
Wukong). 


with  WBL  (k^\  solid  curve)  than  without  WBL 
dashed  curve).  The  relative  decrease  of  the  transfer 
velocity  due  to  WBL  is  calculated  by 


1  = 


EE 

i  j 


EE 

i  j 


(24) 


Rex=u*zo/v,  (23) 

is  the  roughness  Reynolds  number. 

The  horizontally  averaged  DMS  transfer  velocity 
( kw )  over  the  northern  SCS  (Lig.  6)  is  evidently  smaller 


The  parameter  %  generally  decreases  with  the  in¬ 
tensity  of  Typhoon  Wukong  (Lig.  7).  As  the  low  pressure 
of  the  typhoon  reaches  its  minimum  value  of  955  hPa  on 
September  8,  the  value  of  %  has  its  minimum  of  - 13. 1%. 


Fig.  6.  Comparison  of  temporally  varying  horizontal  averaged 
waterside  transfer  velocity  kw  without  WBL  (dotted  curve)  and  with 
WBL  (solid  curve)  for  the  northern  SCS  (15°-25°  N)  during  5-10 
September,  2000. 


8.  Conclusions 

(1)  Effects  of  wave  boundary  layer  on  the  drag  co¬ 
efficient  and  the  sea-to-air  transfer  velocity  of 
dimethylsulfide  (DMS)  are  investigated  using  the 
third  generation  wave  model  (Wavewatch-III)  im¬ 
plemented  in  the  South  China  Sea  during  passage 
of  Typhoon  Wukong  2000.  The  model  was  veri¬ 
fied  using  the  significant  wave  height  measured 
by  TOPEX/POSEIDON  altimeter  for  2000  (Chu 
et  al.,  2004).  Without  the  wave  boundary  layer, 
the  drag  coefficient  and  the  DMS  transfer  velocity 
are  independent  on  the  wave  spectra.  With  the 
wave  boundary  layer,  the  drag  coefficient  and  the 
sea-to-air  DMS  transfer  velocity  are  dependent  on 
the  wave  spectra. 
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(2)  The  drag  coefficient  is  larger  with  the  wave 
boundary  layer  than  without  the  wave  boundary 
layer.  Such  a  drag  coefficient  difference  increases 
with  the  intensity  of  the  typhoon  (in  turns  with  the 
wave  energy).  It  is  only  4%  on  September  5  when 
the  minimum  pressure  of  Typhoon  Wukong  is 
998  hPa,  and  monotonically  increases  to  16.5% 
on  September  8  when  the  minimum  pressure  of 
Typhoon  Wukong  is  955  hPa. 

(3)  The  DMS  transfer  velocity  is  smaller  with  the 
wave  boundary  layer  than  without  the  wave 
boundary  layer.  Such  transfer  velocity  difference 
enhances  with  the  typhoon  intensity  (in  turns  with 
the  wave  energy).  It  is  only  8.3%  on  September  5 
when  the  minimum  pressure  of  Typhoon  Wukong 
is  998  hPa,  and  increases  to  13.1%  on  September 
8  when  the  minimum  pressure  of  Typhoon 
Wukong  is  955  hPa. 

(4)  For  a  modest  typhoon  such  as  Wukong  2000 
(maximum  sustained  winds  of  38  m  s-1),  the 
wave  boundary  layer  shows  the  evident  impacts 
(i.e.,  enhancing  the  momentum  and  reducing  the 
sea-to-air  DMS  flux).  These  features  are  identified 
using  Chalikov’s  (1995)  parameterization  of  the 
wave  boundary  layer  and  McGillis  et  al.’s  (2000) 
parameterization  of  waterside  transfer  velocity. 
The  same  results  (not  shown  here)  are  also  ob¬ 
tained  using  other  parameterizations  (e.g.,  Asman 
et  al.,  1994).  Therefore,  further  investigation  is 
needed  on  the  effects  of  the  wave  boundary  layer 
on  atmospheric  DMS  contents  and  global  climate 
system. 
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